PROCEEDINGS 

°F SCIENCE 



QCD confinement and chiral crossovers, two critical 
points? 



Pedro Bicudo* ^ 

CFTP, Dep. Fisica, Instituto Superior Tecnico, 
Av. Rovisco Pais, 1049-001 Lisboa, Portugal 



E-mail: bicudo@ist.utl.pt 



We study the QCD phase diagram, in particular we study the critical points of the two main QCD 
phase transitions, confinement and chiral symmetry breaking. Confinement drives chiral symme- 
try breaking, and, due to the finite quark mass, at small density both transitions are a crossover, 
while they are a first or second order phase transition in large density. We study the QCD phase di- 
agram with a quark potential model including both confinement and chiral symmetry. We present 
the confinement of quark fields into a flux tube, obtained in SU(3) quenched lattice QCD, to 
iUustrate the importance of the confining quark-antiquark potential. Our finite temperature con- 
fining potential is extracted from the Lattice QCD data of the Bielefeld group. Our formalism, the 
Coulomb gauge hamiltonian of QCD, is presently the only one able to microscopically include 
both a quark-antiquark confining potential and a vacuum condensate of quark-antiquark pairs. 
This model is able to address all the excited hadrons, and chiral symmetry breaking, at the same 
token. Our order parameters are the Polyakov loop and the quark mass gap. We address how the 
quark masses affect the critical point location in the phase diagram. 
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Figure 1: A sketch of the QCD Phase Diagram, according to the collaboration CBM at FAIR [[l|], it is 
usually assumed that the critical point for deconfinement coincides with the critical point for chiral symmetry 
restoration. 

1. Introduction 

Our main motivation is to contribute to understand the QCD phase diagram [p, for finite T and 
/i, to be studied at LHC, RHIC and FAIR, and sketched in Fig[l[ Moreover, our formalism, in the 
Coulomb gauge hamiltonian formalism of QCD, is presently the only one able to microscopically 
include both a quark-antiquark confining potential and a vacuum condensate of quark-antiquark 
pairs. This model is able to address excited hadrons as in Fig. ^ and chiral symmetry at the same 
token, and we recently suggested that the infrared enhancement of the quark mass can be observed 
in the excited baryon spectrum at CBELSA and at JLAB ||]. Thus the present work, not only 
addresses the QCD phase diagram, but it also constitutes the first step to allow us in the future to 
extend the computation of nay hadron spetrum, say the Fig. |T](Ieft) computed in reference [||], to 
finite T . 

Here we address the finite temperature string tension, the quark mass gap for a finite current 
quark mass and temperature, and the deconfinement and chiral restoration crossovers. We conclude 
on the separation of the critical point for chiral symmetry restoration from the critical point for 
deconfinement. 

2. Fits for the finite T string tension from the Lattice QCD energy Fi 

At vanishing temperature T = 0, the confinement, due to the formation of chromo electric 
and chromo magnetic flux tubes as shown in Figs. || and ^ can be modelled by a string tension. 
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Figure 2: First calculation of excited baryons with a chiral invariant quark model. 



dominant at moderate distances, 

^W-l^ + Vo + ar. (2.1) 

At short distances we have the Luscher or Nambu-Gotto Coulomb due to the string vibration + 
the OGE coulomb, however the Coulomb is not important for chiral symmetry breaking. At finite 
temperature the string tension cj(r) should also dominate chiral symmetry breaking, and thus one 
of our crucial steps here is the fit of the string tension o{T) from the Lattice QCD data of the 
Bielefeld Lattice QCD group, [§ |||[ |]. 

The Polyakov loop is a gluonic path, closed in the imaginary time ta, (proportional to the 
inverse temperature T^^) direction in QCD discretized in a periodic boundary euclidian Lattice. It 
measures the free energy F of one or more static quarks, 

P{{))=Ne-^"l'^ , P''(0)P""(r) =A^£'-^'"W/^ . (2.2) 

If we consider a single solitary quark in the universe, in the confining phase, his string will travel as 
far as needed to connect the quark to an antiquark, resulting in an infinite energy F. Thus the 1 quark 
Polyakov loop P is a frequently used order parameter for deconfinement. With the string tension 
o{T) extracted from the qq pair of Polyakov loops we can also estimate the 1 quark Polyakov loop 
P(0). At finite T , we use as thermodynamic potentials the free energy F\ and the internal energy 
Ui , computed in Lattice QCD with the Polyakov loop [glSl, §, 0, |], and illustrated in Fig. |. They 
are related to the static potential V{r) = —fdr with F\{r) = —fdr — SdT adequate for isothermic 
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Figure 3: The fields generated by static hybrid sources, notice how the flux organizes into fundamental 
colour triplet flux tubes. 

transformations. In Fig. ^ we extract the string tensions o{T) from the free energy F\ {T) computed 
by the Bielefeld group, and we also include string tensions previously computed by the Bielefeld 
group [|]. 

We also find an ansatz for the string tension curve, among the order parameter curves of other 
physical systems related to confinement, i. e. in ferromagnetic materials, in the Ising model, in 
superconductors either in the BCS model or in the Ginzburg-Landau model, or in string models, 
to suggest ansatze for the string tension curve. We find that the order parameter curve that best 



fits our string tension curve is the spontaneous magnetization of a ferromagnet []10|], solution of the 
algebraic equation, 

'Tc M 



M 



tanh 



(2.3) 



In Fig. ^ we show the solution of Eq. 2.3 obtained with the fixed point expansion, and compare it 
with the string tensions computed from lattice QCD data. 

3. The mass gap equation with finite T and finite current quark mass mo. 

Now, the critical point occurs when the phase transition changes to a crossover, and the 
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Figure 4: The perpendicular profile of the longitudinal electric field, in units of the lattice spacing a — 
0.07261 (85)/m. The flux tube is so thin, that in quark models it is usually modelled by a single parameter, 
the string tension a. 



crossover in QCD is produced by the finite current quark mass mO, since it affects the order pa- 
rameters P or o, and the mass gap m(0) or the quark condensate (^^7). Moreover, using as order 
parameter the mass gap, i. e. the quark mass at vanishing moment a finite quark mass transforms 
the chiral symmetry breaking from a phase transition into a crossover. For the study of the QCD 
phase diagram it thus is relevant to determine how the current quark mass affects chiral symme- 
try breaking, in particular we study in detail the effect of the finite current quark mass on chiral 
symmetry breaking, in the framework of truncated Coulomb gauge QCD with a linear confining 
quark-antiquark potential. In the chiral limit of massless current quarks, the breaking of chiral 
symmetry is spontaneous. But for a finite current quark mass, some dynamical symmetry breaking 
continues to add to the explicit breaking caused by the quark mass. The mass gap equation at the 
ladder/rainbow truncation of Coulomb Gauge QCD in equal time reads. 



m{p) 
h{p,k,li) 
h{p,k,ii) 



o 

mo + —2 
pi 



dk lA{p,k, lJ.)m{k)p — lB{p,k,n)m{p)k 



271 
pk 



^B+m{kY 
pk 



(3.1) 



(p-k^ + H^ {p + k)^ + H^ 

pk pk \ {p-kf + ^^ 

{p-kY+n^^ {p+kf+^i^^ 2 ^{p+kY+^i^ 



The mass gap equation (Oh for the running mass m{p) is a non-linear integral equation with a nasty 



cancellation of Infrared divergences [pll, 12, 13 1. We devise a new method with a rational ansatz, 



and with relaxation [14|, to get a maximum precision in the IR where the equation extremely large 
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Figure 5: The Bielefeld free Fi energy at T < T^. 




Figure 6: Detail of the string tension fit in the case of T = 0.947^. We cut the low distance part in such a 
way that a linear fit is stable for cuttof changes. 



cancellations occur. Since the current quark masses of the six standard flavours m, d, s, c, b, t span 
over five orders of magnitude from 1.5 MeV to 171 GeV, we develop an accurate numerical method 
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to Study the running quark mass gap and the quark vacuum energy density from very small to very 
large current quark masses. The solution m{p) is shown in Fig. ^ for a vanishing momentum p = 0. 

At finite T, one only has to change the string tension to the finite T string tension o{T) [p^], 
and also to replace an integral in ft) by a sum in Matsubara Frequencies. Both are equivalent to 
a reduction in the string tension, a — )• a* and thus all we have to do is to solve the mass gap 
equation in units of a* . The results are depicted in Fig. ||. Thus at vanishing niQ we have a chiral 
symmetry phase transition, and at finite niQ we have a crossover, that gets weaker and weaker when 
niQ increases. This is also sketched in Fig. |[ 

4. Chiral symmetry and confinement crossovers with a finite current quark mass 

We study whether the two main phase transitions in the QCD phase diagram, confinement and 
chiral symmetry breaking, have two different two critical points or a coincident one. Confinement 
drives chiral symmetry breaking, and at small density both transitions are a crossover, and not a 
first or second order phase transition due to the finite quark mass. 

However the quark mass affects differently these two phase transitions in the QCD. In what 
concerns confinement, the linear confining quark-antiquark potential saturates when the string 
breaks at the threshold for the creation of a quark-antiquark pair. Thus the free energy F{0) of 
a single static quark is not infinite, it is the energy of the string saturation. The saturation energy is 
of the order of the mass of a meson i. e. of 2mo. For the Polyakov loop we get, 

/'(O) -A^e^^"""/^ . (4.1) 

Thus at infinite mq we have a confining phase transition, while at finite mo we have a crossover, 
that gets weaker and weaker when mo decreases. This is sketched in Fig. ^ 

Since the finite current quark mass affects in opposite ways the crossover for confinement and 
the one for chiral symmetry, we conjecture that at finite T and jJ. there are not only one but two 



7 



Two critical points 



Pedro Bicudo 



m(0) 



1 AA 
iUU 


• 




• 


10 


• 




• 


1 1 1 4 


» 1 


10"^ 0.001 0.01 0.1 • 


10 100 


• 








0.1 




0.01 





mO 



m(0) 




mo flnite 



mo=0 



Tc 



Figure 8: (top) The mass gap m(0) solution of as a function of the quark current mass mo, in units of a = 1. 
(bottom) Sketch of the effect of otq on the crossover versus phase transition of chiral restoration at finite T . 



critical points (a point where a crossover separates from a phase transition). Since for the light u 
and d quarks the current mass mo is small, we expect the crossover for chiral symmetry restoration 
critical to be closer to the /i = vertical axis, and the crossover for deconfinement to go deeper 
into the finite pL region of the critical curve in the QCD phase diagram depicted in Fig. [T[ 
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Tc T 

Figure 9: Sketches of the saturation of confinement (top), and of the corresponding crossover in the order 
parameter P polyakov loop (bottom). 

5. Outlook 

We soon plan to complete the part of this work which is only sketched here, i. e. to compute 
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the crossover curves for the chiral symmetry restoration and for the deconfinement at finite T. This 
will also require the study of light hadrons at finite temperature T . Thus we also plan to address the 
excited hadron spectrum at finite T , continuing the work initiated with Tim Van Cauteren, Marco 
Cardoso, Nuno Cardoso and Felipe Llanes-Estrada [^. 
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